Unitarity analysis of general Born-Infeld gravity theories 



ibrahim GiilliiQ Tahsin Cagri JjismanQ and Bayram TekiE0 
Department of Physics, 
Middle East Technical University, 
06531, Ankara, Turkey 
(Dated: December 9, 2010) 

We develop techniques of analyzing the unitarity of general Born-Infeld gravity actions in 
D-dimensional spacetimcs. The dcterminantal form of the action allows us to find a com- 
pact expression quadratic in the metric fluctuations around constant curvature backgrounds. 
This is highly nontrivial since for the Born-Infeld actions, in principle, infinitely many terms 
in the curvature expansion should contribute to the quadratic action in the metric fluctu- 
ations around constant curvature backgrounds, which would render the unitarity analysis 
intractable. Moreover in even dimensions, unitarity of the theory depends only on finite 
number of terms built from the powers of the curvature tensor. We apply our techniques to 
some four-dimensional examples. 
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I. INTRODUCTION 



Tree-level unitarity analysis, that is tachyon and ghost freedom, of a generic gravity model with 
arbitrary powers of the curvature tensors around a constant (nonzero) curvature background is a 
nontrivial problem. On the other hand, for flat backgrounds, only the quadratic terms contribute 
to the propagators, and therefore the analysis is rather simple. In fact, in four dimensions the only 
unitary model, apart from the Einstein's gravity, is the R + aR 2 theory at the quadratic order. 
But, this model is not renormalizable without a fiR\ v term which, when augmented to the action, 
ruins unitarity by introducing a massive ghost mode 

Experience from quantum field theory dictates that at high energies Einstein's gravity should 
be replaced with a theory that has higher powers of various curvature tensors symbolically written 
in the form 

/ = j <fx | i (R - 2A ) + f] a n (Riem, Ric, R, VRiem, . . . ) n | . (1) 

The main nontrivial question is how to find the correct couplings a n that yield a viable unitary 
theory. One might view gravity as a low energy approximation to a microscopic theory such 
as string theory and thus expect to find a unitary (but not necessarily renormalizable) gravity 
theory to any desired order in the curvature by perturbatively computing a n . Of course beyond 
quadratic order this is a very difficult computational problem. Another approach is the so called 
asymptotically safe gravity which conjectures that the dimensionless versions of all the coupling 
constants in ([I]) have a nontrivial UV fixed point and even for infinitely many coupling constants 
the theory has predictive power since the critical surface is finite dimensional 0-0] • In this work, 
encouraged by our recent observation in three dimensions that we briefly summarize below, we 
take a different route and propose that certain Born-Infeld (BI) type gravity actions might define 
unitary models to all orders. Unitarity analysis around constant curvature backgrounds is itself 
a complicated problem when many powers of curvature tensors are involved, here we develop the 
techniques of carrying out this analysis in detail and provide two nonunitary examples in four 
dimensions. In subsequent work [5(, we will give more examples in three and four dimensions that 
are unitary. 

Let us recapitulate the three-dimensional BI gravity action 0] and its success of reproducing 
the known viable theories: 



r 4m2 / r! 
JBINMG = 2~ I d x 



(2) 



where the components of the matrix G read as G^ w = R^ v — ^g^R. In ([2j), we have referred to 
this model as the BI extended new massive gravity (BINMG), since in small curvature expansion, 
this model reproduces the cosmological Einstein-Hilbert theory at the first order, the new massive 
gravity theory 0], which is unitary [§-11], at the second order and the extended new massive 
gravity based on the existence of the holographic c-functions at the cubic and fourth orders 
With the help of the techniques we develop below, we have shown that BINMG is a unitary theory 
at all orders around flat and constant curvature vacua 0] . One of course would like to find analogs 
of (|2|) in higher, especially in four, dimensions. To be able to do this, one has to first establish 
tools for the unitarity analysis which is the purpose of this work. In what follows, for the sake of 
generality, we will keep the discussion in D dimensions and for generic BI actions with the only 
restriction that they reproduce the (cosmological) Einstein-Hilbert theory at the first order. 

The historyof the Bl-type actions is quite rich and for the nongravitational cases a nice review 
was given in [13(]. As for gravity, Bl-type gravitational actions actually precedes a decade their 
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counterparts in electrodynamics. It was Eddington who first proposed that, at least in the absence 
of matter, using the connection as the independent variable, Einstein-Hilbert action can be replaced 
by I = f d 4 x^J 'det R^ v (T) (Note that one actually has to dig this result out from Eddington's 
book, since it is not clearly stated in one place. But, Schrodinger, attributing to Eddington, writes 
this action explicitly on page 113 of his book 11 510 . More recently, Eddington's approach (in fact 
a slight modification of it) was resuscitated in [16] (and the references therein) as an alternative 
to Big Bang cosmology without an initial singularity and with finite density. In 17], instead of 
Eddington's Palatini formulation, the metric formulation where the metric is the only independent 
variable was used in the form I 



f d 4 xydet (g^u + aR^ u + X^ v ) and constraints such as ghost 
freedom on Bl-type gravity actions was studied. Our work follows this line of thought and extends 
the unitarity analysis to constant curvature spaces. We would like to point out to some related 
works where Bl-type gravities, their cosmological and other solutions have been studied [1814221] . 

The main idea of this work is to find a way to obtain the quadratic action in the metric 
perturbation of a generic BI gravity around its constant curvature vacuum, and this can be achieved 
either by explicitly calculating the O (hfa^j action or by finding the equivalent quadratic action 
in the curvature that has the same propagator with the original action. Once the equivalent 
quadratic theory is known unitarity analysis follows with the conventional methods described in 
[l0(. To facilitate understanding and show what is to be expected, let us give one of our results 
here. Let A^ u be an arbitrary (0, 2) tensor built from the curvature tensors, then we will show 

that, in four dimensions, at O {h^ v ) and O (h'^J'j, the action 



i-± 

Ka 



d 4 x 



det (sw + Au, v ) - (aA + 1) \J - det g 



is equivalent to the simpler action 



l O(A>) 



— I d 4 x- 
Ka 



A 



2aA + -A* 



2 ^ 



(3) 



(4) 



where A is the trace of An V . Once this is done unitarity analysis can be carried out with the known 
methods which we shall not repeat in this work. 

The layout of the paper is as follows: In Sec. II, second order expansions of the relevant tensors 
in the metric perturbation hn V are given. Section III is the bulk of the paper which contains 
our general analysis of BI gravities and the corresponding equivalent quadratic actions. We also 
give two examples in four dimensions in this section. Some technical details are delegated to the 
Appendices. 



II. SECOND ORDER EXPANSIONS OF CURVATURE TENSORS 

In order to study the fluctuations of generic BI actions around constant curvature backgrounds, 
we will need to expand various tensors up to second order in the metric perturbation ha V which is 
defined as 

9fj,u = 9^v + Th^ u , (5) 

where we introduced a small (dimensionless) parameter r and a background metric g^ u which 
is quite generic at this stage (i.e. not necessarily constant curvature). [Taking the risk of being 
pedantic, let us note that ([S]) is exact, and that there of course does not exist a natural dimensionless 
parameter in gravity at all scales. So, what one actually means by ([5]) is that in some frame h^ u is 
small compared to g^ for all points in the spacetime, and since there will be another expansion, 
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that is the curvature expansion, r is introduced to keep track of the hn V orders.] Some of the 
computations in this section are actually somewhat tedious but straightforward. They could also 
be found in the literature, albeit somewhat scattered, and probably not in the form we present 
here which proved quite handy in our calculations that follow in the remainder of this work. The 
inverse metric can be found as 

The trace of the metric perturbation is given as h = g^h^. By using these results, the second 
order expansion of the Christoffel connection becomes 

r2* = i%, + r (i*,) l - r% (r^) £ + o (r 3 ) , (7) 

where f p is a background metric compatible connection V p g pv = and the linearized connection 
(r£„) is defined as 

L = \g px {y»K x + v, V - v A v) • (8) 

The main object to consider is the Riemann tensor from which all the other curvature tensors 
and scalars follow. Hence, substitution of = T p „, + 5T P „, to the Riemann tensor R^„ n „ = 

' (JLLS fJLLS fJLLS 1/ fJU 

d p T% u + rJ A r^ - p ^ a yields 

R\ P a = #V + v P (STM - v CT (ai*) + sr^sr^ - sv^sv^, (9) 

where ST^ = r (j^ftJ) L ~ t2 ^$ {^pv) L at trns order. Therefore, the Riemann tensor becomes 

R^pa =R\ P „ + r (R^p*) L - r 2 h» (B? vpa ) L 

- r 2 r a m {(r; a ) L (tQ l - (TZ a ) L (l*)J + O (r 3 ) . (10) 

Note that raising and lowering is done by g pu , but in the above expression, for the sake of notational 
clarity, we do not raise and lower the indices of the linearized Christoffel connection. Here, the 
linearized Riemann tensor (^R IJ, l/pa ^ is defined as 

(#V) L = \ {Vp^aK + v„v„/# - v p v»h au - v a v P K - v a v u w p + v CT v^ v) • ( n ) 

With this result, the second order expansion of the Ricci tensor and the scalar curvature, respec- 
tively, take the following forms 



Rue =Rua + T (R UU ) L - T 2 h^ [R 



V P (7 J 



(r^) L (e) L -(r^(ryj + o(r 3 ), 



(12) 



R =R + T R L + t 2 {R px h ap ht - h v ° {R va ) L - ? a h% 



-va-ua-z 
-9 3 9/3^ 



paJ L 



aaJ L 



r WJ} + °( r3 )' ^ 

where the linearized Ricci tensor and the linearized scalar curvature are defined, respectively, as 

RL = \ i^pS.K + ^iSivK - n Ku - v a v u h) , (14) 



R L = g a/3 Ra(3 - R^Kp- 



(15) 



Note again that the above formulae work for any background space including constant curvature 
spaces which we shall concentrate below. 
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III. BI-TYPE ACTIONS AT O (h 2 ^) 
A. General analysis 

A generic Born-Infeld type action which reproduces the Einstein-Hilbert theory with a bare 
cosmological constant (Ao) at the first order in small curvature expansion is of the form 



I = — I d D x J- det (g^ + A^ v ) - (aA + 1) y/-detg 



(16) 



where A^ u should read as A^ u = a [R^v + f3R^ +0 (i? 2 ) with the definition R^ u = R^ u — j^g^R. 

The O (R 2 ) terms may involve rank (0, 2) combinations of the Riemann and the Ricci tensors, the 
metric and the scalar curvature. It could also involve the derivatives of these tensors, but we 
will not explicitly consider such actions, and we will demand parity invariance, so we do not use 
the a " tensor in the construction of A^ v . Of course, all these technical restrictions can be 
removed and the following discussion can be extended without much difficulty to cover the type of 



actions used in 23[ . Here, the dimensionful parameter a with a (mass) dimension appears only 



beyond the Einstein-Hilbert theory, and k is related to the Newton's constant. Note that in the 



Born-Infeld extension of Maxwell's theory, y — det (g^ + bF^), one must introduce a dimensionful 
parameter 6, since Maxwell's theory is scale invariant, but the BI theory cannot be. On the other 
hand, gravity is not scale invariant and in principle one need not introduce a new scale, one can 
simply use the already existing two scales n and Ao- Nevertheless, introducing a new scale a gives 
more flexibility to the theory. 

To study the unitarity of (|16p . one should consider the quadratic fluctuations around a critical 
point of the action. Assuming that g^ is the critical point and is the fluctuation, we should 
compute the O (h 2 ) terms in the action. To do this by just pulling out the volume density, it is 
convenient to write the action in the form 



2 



I = — I d D x ^g 



det (<5£ + gP^A^) - («A + 1) 



(17) 



Using the second order expansion of the inverse metric, ([6]), and assuming an expansion of A^ u in 
the metric perturbation as 

V = V + ^ff + t* a $ + ( r3 ) > ( 18 ) 

one has 

g^A^ = g^A^ + r (g»>A$ - h^A^) + r 2 (g<*A$ - h^A$ + h^h^ u ) . (19) 

In order to find the second order action in metric perturbation, let us separate the background part 
of gf^A^v and define tB? = g^A^ — g p ^A^ u ^ whose introduction will make the expansion more 
transparent. For a maximally symmetric constant curvature background, one has A^ u = ag^ u where 
a is a dimensionless constant fixed in the theory in terms of the dimensionful parameters such as 
Ao, a, etc. The effective cosmological constant A will also be fixed by the dimensionful parameters. 
For complicated actions, even finding A is a nontrivial problem. The obvious and the conventional 
method is to find the equations of motion and insert the maximally symmetric solution. But, 
finding the equations of motion for these actions is simply too complicated. Therefore, we will give 
a method which bypasses this. Then, becomes 

Bi = (?"A$ - ahi) + r (g^A$ - h^A$ + ah^hjj . (20) 
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Now, we can re-express the BI action with the help of the B£ tensor 
I=JLjd D x | y/- det [(1 + a) 5t + tB p v \ - (qA + 1) j 
d D xy[=g~\ (1 + a) 2 



(1 + a) 



D-4 
2 



det 



(1 + a) 



(1 + a) » (aAo + 1) 



(21) 



where a ^ -1 which is required in order to have a well-defined leading order: if this requirement 
is not put, then the flat space limit cannot be reproduced in the limit of vanishing cosmological 
constant. [For example, if one had fixed a = — with A^ v = aR^, then one would not have a 

proper flat space limit.] Here, the factor (1 + a) 2 is left in front of the determinantal part in order 
not to introduce a factors in the second order terms coming from the expansion of the determinant. 
To find the second order expansion of the action in the metric perturbation, let us Taylor expand 
the determinant in terms of traces up to the order that we shall need 



[det (1 + M)] 1/2 =1 + -TrM + - (TrM) 2 - ^Tr (m 2 ) 



1, 



1. 



1 



6 

With this formula, the second order expansion of \/—g becomes 

1 



48 



+ -Tr ( M 6 ) — -Tr ( M l ) TrM + — (TrM)" 3 + O ( M 4 



(22) 



det g^u 



det (g^ + rh 



-g 



l + T -h + V ( h 



1 K» 



+ T 



(23) 



Then, after using the expansions of the Lagrangian and \J— g in (12 ip one obtains up to O (t 3 ) 



2 D-4 

1= — (1 + a) — 

KOL 



d D x 



(1 + a) 2 -(1 + a)— (qAq + 1) 



T 



(l + a)BP+ (l + a) 2 -(l + a)— (aAo + 1) 



r-2 r 



BP) 2 - 2BP v B v p + 2 (1 + a) hB? + [(1 + a) 2 - (1 + a)^ (aA + 1)] (h 2 - 2a 2 ,) 



(24) 



O (t ) term just gives the value of the action for the vacuum solution and it will not be relevant 
anymore. But, it gives us some crucial information about the Bl-type actions, that is for even 
dimensions the value of the constant curvature is not bounded by the action; however, for odd 
dimensions a > — 1 is required for the reality of the action. Now, we would like to go back to our 
original tensor A^ v . First, we write the O (r) term in the above expression in terms of A^ v . This 
term gives the nonlinear equation of motion for the constant curvature background. One needs to 
find what the zeroth order of B p p is in terms of A^. This is given as 



BP = g^A^ p -ah + 0{T). 



Then, the action at the first order reads 



l O(h) 



(1 + a)" 



d D x 



(1 + a) (g p> *A$ + h) - (1 + a)^ 1 (aA + 1) h 



(25) 



(26) 
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After removing possible boundary terms, taking the variation with respect to h^ u or more concisely 
looking at the coefficient of h^ u and equating it to zero yields the source-free nonlinear equation 
of motion for a constant curvature background, namely, the equation of motion that relates A to 
Ao and the other parameters of the theory. Hence, to get the vacuum of the theory, one need not 

explicitly find the equations of motion which is straightforward but quite tedious. 

2 



Now, let us find the quadratic action in h pv in terms of A„„ 



terms in 



can be written in terms of A^ v as 



.The (bp) -2BP B» + 2(1 + a) hBP 



(BP)" - 2B£B" p + 2 (1 + a) hBP = (g^A$ 



2A^ A^ v 



+ h^ 



AaA$ + 2g^gP°A$ 



2a 2 h 



[IV 



a(2 + a)g^h . (27) 



Contribution coming from the r (1 + a) BP term in (I24p is 

BP p = O (r°) + r [g»»A$ - (A$ - ah, u ) 
In all together, the quadratic action in h^ v in terms of A pv boils down to 

(1 + ° )V f d °^{ \A$A^ - \ (g^AWf - (1 + a) fA® 



(28) 



l O(h 2 ) 



KOL 



1 



1 r 



l-(l + a)^(aAo + l)] (h 2 -2hl„)X. (29) 



+hT {A® - -9,u9 pa A p 1 J J ) 

To remove a possible confusion coming from the notation, we should note what is represented by 
the term A^A^y. It is basically A$A^ = g^g v ^ A^u A^, that is A^ does not represent the 



first order of A^ u . If required, we show the first order of A^ v as (A^)^. Equation (|29p is our 
main formula which can be applied to any Bl-type action for any value of the constant curvature 
[i.e. we have not done a small curvature expansion, that is, the formula at O (h 2 ) takes care of 
all the contributions coming from all powers of the curvature] . Let us summarize what one needs 



to do to analyze the unitarity of a given BI gravity: One computes A$ and A£J, and using (f26l) 
one finds the vacuum of the theory, and finally computes the O (h 2 ) action via (|29p . Then, this 
action can be studied using conventional techniques that were discussed in [lOj . Of course, as we 
shall see below with some examples, depending on the complexity of A^, explicit computation 
of (I29p could be a very cumbersome problem in generic dimensions. But, a close scrutiny of 
it reveals remarkable simplifications in even dimensions, higher than two, and especially in four 
dimensions. Such simplifications, in four dimensions, will provide us with another method of 
analyzing the unitarity of the BI gravities, namely, the method of Hindawi et al 24] that leads to 
the construction of an equivalent quadratic action (in curvature) whose unitarity has been already 
studied by conventional methods. Let us concentrate on D = 4 first whose action is 



i (2) 



l O(h 2 ) 



1 



-(J 



4^(i) " I [sTA^f - (1 + a) ?*A® 



( Aft) - - 



fx* AW 
per 



1 . 

+ |OA 



h z - 2h 



fW 



(30) 



By examining this action, one can figure out an interesting relation between the metric pertur- 
bation expansion that led to this action and the A pu expansion of (|16j) . Remember that A pv is 
dimensionless, so assuming proper convergence, a Taylor series expansion over A pv is legitimate. 
If An, v involves terms of O (R 2 ) and/or any other higher curvature terms, the A pu expansion is 
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not simply equal to the curvature expansion in which the expansion is over the nondimensional 
quantity aR. Let us write symbolically the expansion of f 1 1 6 [) in A^ v as 



I = — I d x 



~ — I d x J—g 

Ka 



det (g^ + A^) - (aA + 1) \/-detg 

d x yj — g 



c n A n - (aA + 1) 

n=0 



1 2 °° 

-(R- 2A ) + — E 
k Ka ~, 

n=2 



, (31) 



where the last equality follows from our assumption that Einstein-Hilbert action is reproduced at 
the lowest order. Note that, up to n = 3, this expansion can be obtained with help of (|22p . and 
the n th order term represented with A n involves terms like A n , A n ~ 2 A 2 lu , A n ~ 3 A^APA'^ l , etc. In 
principle, each order in (|3ip contributes to the quadratic action in the metric perturbation given in 
(|30|) . but we will see that this is not the case in four dimensions. The O {h 2 ) contributions coming 
from the O (A n ) terms where n > 2 have the form 



r(n) 
'o(/j 2 ) 



d 4 x 



gc n - 



A 



n-2 



d n iA$A^ + d n2 (g^AW 



in— 1 



d n3 A$ + (Uh^AM + d n5 hg pa A$] + A n 



pa 



dnsh^v + d n -jh 2 



(32) 



where A is defined as A 
bers. Therefore, the O (h 



fiu = ag^u as above, and the coefficients d n are just num- 
2 ^ contributions coming from the O (A n ) terms are in the form of 



[e n 2 (h) a n 2 + e n i (h) a n 1 + e„o (h) a n ] . Hence, one expects that if each order in the A^ u ex- 



pansion of (I16p contributes to the quadratic action in metric fluctuations, then that action will be 
composed of the seven terms specified in (|32p with a coefficient which is a power series in a. With 



this result, one can trace the contribution coming from each order in (|31f) to the O (h) action 
(|3Up . Let us investigate each term in (|30p in order to find which orders in the A^ expansion 
contributes. The first two terms in (|30p . which are quadratic in A^, have coefficients that do not 
depend on a. Therefore, these two terms involve O (h 2 ) contributions only coming from the second 
order terms in the An V expansion of (|16|) . The coefficient of the third term in (|30|) is (1 + a), so 



it is composed of contributions coming from O (A) and O (A ) terms in the A^ expansion (131 p . 
The fourth term has a coefficient which does not depend on a, so it comes from the first order 
of the A^ v expansion. Thus, all the O (h 2 ) contributions coming from O (A n ) terms with n > 2 
are identically zero for four-dimensional Bl-type actions, and as we will see this curious case has a 
generalization to higher even dimensions. With these observations, one can deduce the fact that 
in four dimensions (|3U|) can be obtained first by making an expansion in A^ u up to third order via 
f|22p . and then by finding the quadratic action in metric fluctuations. In other words, remarkably 
the free theory of the following actions are exactly the same: 



2 

Ka 



d A x 



det (g^ + A M „) - (qA + 1) \J - det g 



and 



o(A 2 ) 



2 

Ka 



d x yf—g 



l - g ^ v " « A o + \ (jgTA^f - \a 2 



[IV 



(33) 



(34) 



which was obtained by expanding (|33p . Here, note that we truncated the A^ expansion at the 
second order, but we do not require A^ u to be small. This truncation can be done and the equality 
of the above two actions at the free level can be achieved merely due to the fact that contributions 



9 



of the higher order terms in the expansion to the quadratic action in metric fluctuations are 
identically zero. Such a remarkable cancellation in four dimensions is related to the fact that we 
have the square root of the determinant of a linear combination of matrix functions one of which is 
expanded around a constant curvature space and it would not work for a generic background. Let 
us verify this result by explicitly calculating the quadratic action in metric fluctuations for (|34p . 
However, to be as general as possible and to see some cancellations, let us work in D dimensions 
where only measure in (|34f) changes to d D x. Then, expanding each term in (|34f) by using (I23p and 
(|19|) with A^y = ag^u one has 



ST V = aD + T (sTA® ~ ah) + r 2 (f»A<$ - hTA$ + ah 
and all together up to quadratic order, the action reads 



(35) 



l O(A 2 ) 



1 



— I d D x 



aD 



+ T 



a 2 D a 2 D 2 
aD 



aAn 



9 / y 9 \ A 



/i — aAnh 



l + f-a) 9 ^A$ 



(36) 



2 

(D-4) 



a + 



D-6 



„>)( 



h 2 -2h^) + 



aK 



°(h 2 



2h 



In obtaining this result, one should rewrite (A IJ,V )^ and (A^)^^ coming from A 2 ^ in (|34p in terms 



of A$ and A$ as 



'(i) 



(AH(i) = (gT^Ajf* = r a 9 vfi A% - 2ahr. 
Let us now concentrate only on the O (r 2 ) terms: 



(38) 



0(A>) 
( 0(/i 2 ) 



1 



1 + r ! 



+ (l + ^ - 2a) ftT (A« - \ 9 ^A^ 



(39) 



(D-4) 



a + 



(£> - 6) a 2 



(tf-2^) + ^(tf-2/&,) 



In four dimensions, 



reduces to (|3U|) as it was promised. In Appendix A, we give a simple 



example with two-dimensional matrix functions that shows the connection between the A^ v ex- 
pansion and the metric perturbation expansion. In generic even dimensions with D = 2n + 2, if 
one wants to carry out a similar analysis, then one has to expand up to O (A n+1 ) with n > 1. But, 
again we should stress that the compact formula (|29|) works all the time without recourse to such 



an expansion. However, depending on the complexity of A^, one can choose to use either the ex- 
pansion method or the compact expression. As for odd dimensions, because of the nonpolynomial 
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D-4 

prefactor (1 + a) 2 in qju uu C .dino ±±± .uc 
contribute. The most efficient way to get the quadratic fluctuations for odd dimensions is to use 



all the terms in the A„ v expansion (or the small curvature expansion) 



A similar analysis can be done for the 0(h) action in four dimension which is 



iO(h) 




(l + a)g^A^ + (a-aA )h 



(40) 



This action involves contributions coming only from the second order expansion of (|16p in A^ v 
just as the O (h 2 ) action. In order to understand this behavior, let us first look at the O (h) 
contributions coming from the O (A n ) term for n > 2 



An) 
2 0(h) 



d 4 x 



gc n 



A" 



~ l d n \ (l 



g^AW)+d n2 A n h 



(41) 



where d n coefficients are just numbers. Therefore, the O (h) contributions coming from the O (A n ) 
terms are simply in the form of [e n i (h) a n_1 + e n o (h) a n ] , since A ~ a. Hence, one expects that 
if each order in the A^ v expansion of (|16|) contributes to the linear action in h^, then it will be 
composed of the two terms specified in 041 p with coefficients that are of the form a n . With this 
result, one can trace the contribution coming from each order in (131 p to the O (h) action (|40p . Let 
us investigate each term in (|40p in order to find which orders in the A^ expansion contribute. 
The first term in (|40p has a coefficient of (1 + a). Therefore, this term involves O (h) contributions 
coming from the second order terms in the A^ v expansion of (|16p . The coefficient of h in (|40p is 
also first order in a, but this time it implies that only the first order of A^ u expansion contributes. 
Therefore, the vacuum of f 1 1 6 [) and (|34|) are the same. One can verify this result explicitly from 
O (h) of the O (AjJ\ action which can be read from (1361) as 



G(A 2 
l O(h) 




1 + 



aD 



a 9 



a(D-2) 



1 + 



{D-4) a 



h — aAoh 



This action reduces to (|40p in four dimensions. Just like the analysis of O (h 2 ), for generic even 
dimensions D = 2n + 2 one has to expand (I16p to O (A n+1 ) with n > 1, then find the vacuum of 
the theory. For odd dimensions, since all the powers of A n contribute, the most efficient way to 
find the vacuum of the theory is to use (|26p . 



B. An example 



To apply our tools, for the sake of simplicity, let us consider the following model which we know 
to be nonunitary even around the flat space: 



d 4 x 



det (g^ + aR^ u ) - (aA + 1) \/-detg 



(42) 



Here, according to our results above, one expects (which we shall verify below with several different 
techniques) that the second order action in the metric perturbation h^ u involves contributions only 
coming from the O (aR) 2 



expansion: 



l O(R 2 ) 



2 

Ka 



d 4 Xy 



-g 



n 



(R - 2A C 



a 

T 



R 



1 



-R< 



(43) 



Therefore, the O 


(aR) 3 ' 


, o 




(h 2 ) (0D and i 


33D are equivalent 



theory of (j4^|) around the extremum of it by using ([3H|) . 
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1. Analyzing the BI action formed by the Ricci tensor via second order perturbations in h 



Let us define Au V = aRu V . Then, A 



terms of Aq. Then, A$ is given as 



>-\AV 



aAg 



a = aA , where A will be determined in 



A (l) = aR L 

(2) (2) 

and Aiip = aRuJ and referring the details to Appendix B, we have 



(44) 



(45) 



First of all, let us determine the nonlinear equations of motion for the constant curvature back- 
ground which will relate A to Aq by using 



hih) = — I d 4 x 



1 



d 4 x a/— 5 



1 + a) (s"M$ + h)- (aA + 1) h 
* (1 + aA) (VuhT - V^/i) + a (A - A ) h 



(46) 



Note that this first order correction should be zero around the extremum, therefore, after dropping 
the first term which is a boundary term, one has A = Aq. As for the second order action, one has 



l O(h 2 ) 



- (a + a 2 A) hT Qi?^ - ^uRl ~ j9^ 



(47) 



a 



where f d 4 x a/ —gR^ u R^ is calculated in Appendix B as 



cTx yZ-gR^R 



f g^D - V^V U + A^J R L + - —g^R L 



2A 



14A r A 



3 ' :'» 

Then, after some algebra the quadratic action reduces to 



8A 2 



(48) 



'a 2a 2 A x 
2 + — 



2A_ 



(A-A ) (V-2^ 



(49) 



where we kept the background gauge noninvariant term (the last part) just to show an intermediate 
step of the computation. Once A = Aq is used, one ends up with 



o(h 2 ) 



1 

an 



a 2q 2 Aq \ a 2 



2 + 3 



7 flU 



°0. 



2A n 



-9/iiyRL 



(50) 
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This action is exactly equivalent to the linearized action one obtains from the O \ (aR) 2 action 
[Note that the linearized version of (j4"3j) has been worked out in several places [10|, |25[ , and we also 
reproduce it below.] The fact that (|50p has at most a 2 terms show that the contributions coming 



from all O 



(aR) 



n+2 



vanish. We stress once again that this is a highly nontrivial cancellation 
brought by the determinantal structure of the action. It is worth to study explicitly how this 
cancellation takes place at O (aR) 3 



2 

not 



1 0{R?) 
and defining 

one has 



a 



which we do now. At this order the action reads 

.,2 / 1 \ „,3 



(R - 2A ) - ^ (r%, - l -R 2 \ + ^ (aRwRPRT - QR^R + R 



K = R% 



-R 2 



S = 8R pu Rp a R a u 



6RRp V + i? 3 , 



0(R 3 ) 



d £a/— g 



(R-2A )-^K + ^S 



(51) 
(52) 

(53) 



Finding the O (h 2 ) action of this theory is a very cumbersome problem. To somewhat simplify 
this, one can first find the equations of motion then linearize the equations of motion and then do 
the reverse calculus of variations procedure to get the action. Of course in this process boundary 
terms are dropped and one has to be careful with an overall sign that can be fixed by coupling the 
gravity action to matter. The equations of motion follow as 



KCt 



a 
4 

a 2 



(R - 2A C 



+ 




T 




a 3 


+ 


4 




a 3 


+ 


8 




a 3 



RR 



2R\vauR Xa — n ( R 



a 

9/J.u + 77 



[IV 



WUfi 



:9^uR 



+ (2R?R pa R« + [ 9lxv V a V/3 (RPPRf) + □ (R P V R W ) - 2V a V M (r p u R«) } ) (54) 



([2V a V M (RR«) - 9^V a Vp (RR a ?) - □ {RR t 

1 



[{g^n-VuV^+R^] (R 2 a p 



-R< 



2 (5 -/matte 



where we defined the energy-momentum tensor as — — -j= 
ground (R^ u = Ag^ u ) should satisfy source- free equations of motion with the results 

1 - 



1RR p v R pp 



Constant curvature back- 



K = R 



-R z 



S — 8R RuaRu 



QRR 2 ^ + R 3 



0. 



(55) 



Then, the equations are satisfied if A = A$. Now, let us linearize (|54p around its vacuum (defining 
Tfiu (h) = 5 (p^-)) by use of the formulae in Appendix C and 



5K 



-2AR L , 



5S = 0. 



(56) 



The linearized equations of motion after using the source-free equation of motion for constant 
curvature background becomes 



T„ u (h) 



1 4aA 
K 3k 



a 



2 k 



2A °- R 1 



(57) 
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which exactly matches the equations that result from the matter coupled version of the action 
as promised. This shows explicitly that O (ai?) 3 terms cancel each other. This cancellation will 
work for any arbitrary order beyond this, as we will show with a different method below. 



2. Another method for unitarity analysis 



Hindawi et al j^] gave another method of analyzing a generic higher derivative gravity model 
by reducing it to the equivalent quadratic curvature theory in the sense that it has the same free 
Lagrangian as the original higher derivative theory. Here, we will review their approach and apply 
it to our example (142 p . Before we describe their method, we should note that unlike our method 
which led to the compact formula (|29p that works in all cases, the method of Hindawi et al works 
only when one deals with not matrices but scalar objects or one has a finite number of curvature 
terms. Keeping this caveat in mind, which will be better understood below, when Hindawi et al 
method works, it provides a fast algorithm in getting the equivalent quadratic action. 

To understand the essence of the Hindawi et al method let us consider the following simpli- 
fied problem. Suppose we have a function f(x(t)), and we would like to find the e 2 order of 
/ (x (to + e)). But, instead of doing this, we can find a function g (x (t)) = ao + a\x (t) + a2X 2 (t) 
whose second order expansion around to yields the same second order expansion of / (x (t)) around 
the same point. After some straightforward analysis, one can show that g(x(t)) can be ob- 
tained by expanding f(x(t)) around xq = x (to) up to and including O (x (t) — xq) 2 , since 



O 



(x (t) - x 



,2+n 



gives e 



corrections with n > 1. Hence, one can read the coefficients for 



the correct g (x (t)) to be 



ao = f(x )- 



dx 



xo+'z 



1 



dx 2 



' 0' 



«1 



X{) 



\df- 




\d 2 f] 


_dx_ 


-Xo 

XQ 


dx 2 



a 2 



dx 2 



(58) 



X() 



Note that if one just wants the "equations of motion," then one carries out the above procedure at 
0(e). In this example, / represents the Lagrangian, x any curvature tensor or scalar, and e rep- 
resents the metric perturbation hnv ( x )- Similarly, to, xq are used in analogy with the background 
metric g^ u , etc. 

a. Cubic theory Now, let us turn to our example (|42|) and to specifically its third order 
expansion in curvature given in (|5ip . In order not to introduce the metric or its inverse during the 
expansion around (^R, , let us take the Lagrangian density of (|5ip to be a function of R and 
R£ as 



a 



f(R,R<i) = -(R-2A Q ) 



a 

T 



RvRp 



l 



r2 ) + is ( 8R p rp »K - SKKR + R 3 ) ■ (59) 



Expanding / (R, R£) around (R, R^J with the assumption of small fluctuations about the back- 
ground yields 

" df 



S'quad-equal 

(R,R%) =f(R,RZ) 
~d 2 f 



+ 



df 



dR 



(rM) 



dR 2 
d 2 f 

dR&dR^ ^ R m) 



(r-r) 

(rM) v ' 



(r-r) 2 + 



+ 



dR- 



dRdR^. 



(RM) 
(R1 - Rf) (rp - K 



\Kp - tip 

(RM) 

(R-R)(R%-R$) (60) 
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Computing the relevant derivatives one ends up with 



<?quad-equal 



a a 2 - a 3 
2 + T R+ ~6 



R 2 - 2RX 



R-R 



+ 



2 

2 

1 fa 2 



+ ^-r] (r-r) 2 



(61) 



■ ^rPUr-r) (r<$-R 



1 

+ 2 



~ 5 X + t { 25 pK + 2RX - sfcR) (r$ - Rf) (K - K 



2 ? a 4 

For constant curvature backgrounds, the corresponding quadratic action becomes 



d x^J —g 



a 



" 2 ^2 a 2 D 2 



-(R-2A ) + —R Z -—R^ 



(62) 



which once again shows that the cubic term in (|59|) does not contribute to the free theory. We 
should stress that if one takes arbitrary coefficients instead of the ones we have which are (8, —6, 1) 
at the cubic order (|59p . then one would get a different quadratic action that does not follow from 
the A^y (in this case it is just aR^) expansion of (j4*2j) . 

Now, let us also obtain the source-free nonlinear equations of motion for a constant curvature 
background by finding the equivalent action at O (R). Similar steps lead to 



5lin-equal 



+ 



Of 



dR 



{rM) 



(r-r) 



+ 



and to the action 



(1 + qA) 



R-2 



A + oA 2 " 
1 + aA 



df 



+ R 



(r1 - Rf) , (63) 



(64) 



where we used R^ = Ag^ u . Then, identifying A = A °^"^ , one obtains A = Aq. 

b. Full nonlinear action We mentioned above that Hindawi et al method does not work when 
one deals directly with matrices. Let us show this with 

/ {Rfj,u 

is anal 

v / detM 



det (g^ u + aApv), (65) 
and try to find df which is needed for this analysis. Defining = g^ v + aA^ u , one has 

M -1 dAfl , 



df = d (V det M 



Tr- 



im) 



where M _1 dM is an ordinary matrix multiplication. Here, the basic problem is to find M _1 which 
cannot be done in exact form for a general A^ u and even when A^ v = R^ u . But, one can always 
expand the determinant in terms of traces and apply the Hindawi et al method. Even though this 
is the case, for a complicated A^ v the determinant will yield many terms in generic dimensions 
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and as we show below even for four dimensions. Let us consider the action 
formula 



and use the exact 



detM = l{( 



(TrM) 4 - 6Tr (m 2 ) (TrM) 2 + 3 [Tr (m 2 )] + 8Tr (m 3 ) TrM - 6Tr (m 1 
for M = 5^ + fti?^ , one gets 



(67) 



2 3 4 

det (<5£ + ai?£) =1 + ai? + y # 2 + y i? 3 + ^i? 4 



^-RR^Rp 



a 



—RliRPRi 
3 p v P 



4 4 4 

a p2 pu pi' i a p D(i pp pf I tt pit p^ pp per ' pu pp per tjv 

— K RyR^ + —KJv^tVyKp + —K^K^K^Kp - — K^K^K u Kp. 



(68) 



Defining / (R, R£) = y det (5% + aRu) and with the help of (|6U|) and the formulae in Appendix D, 
we have the corresponding quadratic action as 



7 



2 



d 4 X\ 



-g 



a: 



a 



(R - 2A ) + —R 



ft 



-i? 



(69) 



Once again we have proven that the O (R 2+n ) with n > 1 terms do not contribute to the free 
theory for the exact BI action (|42l) around its constant curvature vacuum. Note that with the help 
of an equivalent action at the linear level as we have done before, 



I 



g \ -(1 + oA) 



R-2 



' A + aA 2 ' 
1 + ftA 



(70) 



setting A = ^f^A > one nas ^ = • We should stress that to get this result with the conventional 
method of finding the field equations and looking for a solution of the form Rp U = Agp U is highly 
cumbersome for an action which is given as the square root of det (5% + aR^) 



C. Unitarity of the theory proposed by Deser and Gibbons 

While constructing the Bl-type gravity actions, among various criteria, one of the easiest to 
realize is the unitarity of the model around flat space. This means when small curvature expan- 
sion is carried out at the quadratic order in four dimensions, one should get the unique theory 
i (R - 2A ) + aR 2 + 7 (i? 2 



4i? 2 j, + Rp\ vpa J which is free of ghosts. Deser and Gibbons 17] sug- 
gested that at the quadratic order, one should get, dropping the aR 2 term, only the Einstein plus 
the Gauss-Bonnet combination (the 7 term). We will study such actions in a separate work, but 
here let us consider an example (the one suggested by Deser and Gibbons) of these models in which 
one does not have quadratic terms when expanded around small curvature: 



d 4 x 



det 



ft / 1 
gp U + aRp U + — [RppRi — -RRp U 



1 



(oA + 1) v 7 - det g 



(71) 



It is easy to see that the lowest order correction to the Einstein-Hilbert theory goes like O (R?) , 
which means around flat space the graviton propagator is the same as that of Einstein-Hilbert 
theory. (Note that for flat space to be the vacuum, one also sets Ao = 0.) But, around its 
constant curvature vacuum unitarity of this model has not been checked before, since it is a highly 
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nontrivial computation without the tools we have developed above. To carry out the analysis, 
we can find the O (A^J^ action which has the same O (/i 2 ) action as (j71j) . Here, A^ u = aR^ u + 

(^R^pRP — ^RR^jA and let us stress again that O (^4jL) action is not equivalent to O (aR) 4 
action. If one naively does the latter expansion, one will simply get an inconclusive result since 
one would have neglected the O (aR) 4+n corrections. But, cin expansion in A-^y tcikcs care of all 
the relevant terms and cancellations. Therefore, using (|22p we have 



l O(A 2 ) 



— I d A x 



-(R-2A ) + — 



RR^R^ 



R^R P R U 

p V fJi 



4 R 



+ 22 [ R v R, il R a R p - t; r2r u r u + -R A 



1 

T 



2R^RPR v a R^ + 2RR^RPR^ 



(72) 



Now, let us just concentrate on the higher curvature terms and define 
v3 / x 



a 



RR^R^ 



R^R P R U 

p V f-l 



a 



+ 32 i R v RV h R o R °p 



-R 2 R p R v il + -R 4 - 2R p R p ,R u a R° l + 2RR^R 9 V R V ^ 



(73) 



The first thing we should find is the correct A which can be found by using the first order expansion 
of / (R, ffl£) around ( R, R 9 ) = (4A, A<5^). This procedure leads to the equivalent linear action 



i-± 

not 



3a 3 A 2 ~ 4 

: + 



crA 3 



R-2 



a 



3 A 3 + 



3a 4 A 4 



(f 



from which one can get the equation that determines A 



An + A + 



o 



2 A 3 a 3 A 4 



3a 3 A 2 , a 4 A 3 

A < A 



0, 



+ 



o(r 2 ) 



(74) 



(75) 



which has real roots, but they are not particularly illuminating to display here. (One thing we 
can note is that even for An = 0, there are two real roots one of which is nonzero with a value 
A ~ 2.59/a.) Now, we can employ the Hindawi et al method to get the equivalent quadratic action 
using (|6U|) and the relevant results of Appendix D: 



KOL 



d^X 




a 



l A 3' 



R 



(76) 



For generic a, this theory is plagued with a massive ghost Thus, the action proposed 

by Deser and Gibbons 17J does not yield a unitary spin-2 theory around its constant curvature 
background for any choice of the curvature except the flat space. But, setting a = 4 one can get 
rid of the "bad" i? 2 ^ term, and hope to obtain a unitary theory. However, this turns out to be not 
true, since in this case setting A = An which follows from (|75p . one ends up with 



/ 



d X\ 



1 



(R - 2A C 



(77) 
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which has the opposite sign of the Einstein-Hilbert action. That means as long as one has k > 
(which we must have for the unitarity in fiat space), the small fluctuations will have negative 
kinetic energy and even for the tuned value of a, (|7ip defines a nonunitary theory. We should note 
in passing that this result does not necessarily imply negative energy for the exact nonvacuum 
solutions such as black holes of (|7ip . We have not yet found the black hole solutions of this action, 
but we can give an example in which small fluctuations around the vacuum have negative energy 
yet the exact solutions have positive energy. This example is the Einstein-Gauss-Bonnet theory 
whose exact spherically symmetric solution was given in [26] and whose energy was computed in 
{25]]. As discussed in the latter work, this energy is positive, even though the linearized action 
of the Einstein-Gauss Bonnet theory around its constant curvature vacuum is opposite to that of 
Einstein's theory [just like fj77j) ] . This is because the spherically symmetric Schwarzschild-de-Sitter 
solution goes (say in five dimensions) as —goo = g rr ~ 1 + ™ + Ar 2 unlike the usual Schwarzschild 
solution which goes like —500 = 9 rr ~ 1 — ?Zj the two minus signs take care of each other. 



IV. CONCLUSION 



We have developed techniques of analyzing the unitarity of Born-Infeld gravity actions around 
their constant curvature vacua. The special determinantal form of the action gave rise to remark- 
able simplifications that allow one to write a compact expression for the free, that is O (h 2 ) , theory. 
To summarize our result, let us note the following: One needs to find the O (h 2 ) action of 



det (5H + aRii + /3 (Riem, Ric, R, ...)£) - (aA + 1) 



(78) 



to study its tree-level unitarity. In this work what we have done is to give a method to determine 
the parameters K, A, a, b, c in the following Lagrangian whose O (h?) expansion equals that of 

(USD 

^equivalent =T77 \ i R ~ 2A ( K ' Z 3 ' A , • • • )] 

K (K,a,p, A , • • • ) 

+ a (k, a, 13, A , . . . ) R 2 + b («, a, /3, A , . . . ) R 2 ^ + c (k, a, /?, A , . . . ) -RjL p(7 , 
to all orders in the curvature expansion. Once this equivalent quadratic Lagrangian is obtained, 



unitarity analysis proceeds with the standard methods as discussed in [lOj]. We have also presented 
two examples one of which was proposed as a unitary theory in flat space 17J , but turned out to 
be nonunitary in curved space according to our computation above. The other simpler example 
was considered to show the details of our method. 

Let us give a recipe of how one should check the tree-level unitarity of a given Born-Infeld 
gravity in generic dimension D around its constant curvature vacuum. First to find the effective 
cosmological constant A, one has to expand the action up to O (h) around the constant curvature 
vacuum using (|26|) . or one should find the equivalent linear, that is O(R), action and read the 
cosmological constant from it. Then, one should find the O (h 2 ) action using ()29j) or alternatively 
one should construct the equivalent quadratic action [O (R 2 )]- The method we have presented (|29|) 
works just as good in odd and even dimensions. But, the second method, as discussed in detail 
in the text, which proceeds by construction of an equivalent O (R) and O (R 2 ) actions should be 
done with great care depending on the number of dimensions and on the complexity of A^. The 
original BI action cannot simply be expanded in small curvature to get these equivalent actions via 
(|22|) . What always works, in principle, is that the determinant can be expanded exactly in terms of 
traces within the square root, then one can do the expansions f!63[) and (|60|) . and use the Hindawi 
et al technique. But, the exact expansion of the determinant in terms of traces can generate quite 
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a large number of terms especially for D > 4. [For example, in (|71|) doing such an exact expansion 
is not advised to the reader.] Therefore, to get the equivalent action one should proceed as follows 
in generic even dimensions D = 2n + 2: one has to expand the BI action up to O (A n+1 ) with 
n > 1 using (|22p . if the resultant action is not already quadratic in the curvature, then using 
([GO)) , the equivalent quadratic action should be constructed. For generic odd dimensions, the best 
way is to use (I29p . but for D = 3 and for not so complicated A^ u , exact trace expansion can also 
be employed. In this work we have laid out the details of checking unitarity of BI gravities, in 
a separate work we will provide examples of unitary models around flat and constant curvature 
backgrounds 0]. 
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Appendix A: A Two-Dimensional Example 



In order to understand why in even dimensions finite number of terms in the expansion of the 
Bl-type actions contribute to O (h) and O (h 2 ) expansions, let us study a simple two-dimensional 
determinantal function 



l{r,l) 



\ 



det 



1 a(r) 6(r) 

1 +7 c(r) d(r) 



(79) 



where r and 7 are two independent variables. The r, 7 expansions of / (r, 7) represent the metric 
perturbation expansion and the A^ u expansion, respectively, for the Bl-type actions. What we will 
show in this Appendix is that / (r, 7) and the function g (r, 7) defined as 



5 (r, 7 ) = l + _ 7 [ (T) + d(r)] 



have the same the O (r) expansion around r = only if 



a(r) b(r) 
c(r) d(r) 



T=0 




(80) 



(81) 



which is the analog of the maximally symmetric constant curvature background in the Bl-type 
gravity. Here, the important point about g (r, 7) is that it is just the O (7) expansion of / (r, 7) 
obtained by using (|22[) . but note that we exactly define g (r, 7) in this way and do not assume that 
7 is small. Thus, staying at first order in r expansion requires just the first order in 7, while one 
naively expects that first order in r expansion should involve each order in 7. Let us understand 
this in more detail by considering a generic function 4> (r, 7) and expand it in r as a Taylor series 
around r = 



^(r, 7 ) = 0(r = O,7)+(^) 



r + O r 2 +..., 



(82) 



where (ff ) is a function of 7 only. One can write the power series expansion of (757) i n 7 
by assuming 4> (t, 7) = X^o "4>i ( r ) 7* anc ^ expanding each (r) to the first order in r. Then, one 
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has 



#(r,7)=X>< (r = 0)7* + 



i=0 



E 

Li=0 



7 



T=0 



rH 



r=0 



E 



0r 



T=0 



y. (83) 



For the determinantal function / (r, 7), the terms (jjjtf^j , i > 2 are all zero. Let us observe this 
for the i = 2 term explicitly. First, one can have the O (j 2 ) expansion of / (r, 7) by using 



\ 



det 



1 
1 



+ 7 



o(r) b(r) 
c(r) d(r) 



l + ^7[«(r)+d(r)] + i7 2 [a(r)+d(r)] 2 



Assuming 



one has 



o(r) 6(r) 
c(r) d(r) 



1 2 

4 7 



a 
a 



a 2 (r) + 26 (r) c (r) + d 2 (r)] + O ( 7 3 ) . (84) 



o(r 2 ), 



(85) 



/ ( r , 7) = 1 + -7 [(00 + roi) + (a + rdi)] + -7 [(a + roi) + (a + rdi)]' 
(a + rai) 2 + 2 (r&i) (rci) + (a + rdi) 2 l + O (7 



"4 7 

(1 + 7 a ) + i 7 r (ai + di) + ^ 7 2 Ua 2 + 4ra (a x + d x ) + O (r 2 



(86) 



V 



2a 2 , + 2ra (a x + di) + O (r 2 )] + O (7 s ) 



= (1 + 7 a ) + -7T (ai + di) + O (r 2 ) + O (7 s 

Thus, O (r) contributions coming from the two O (7 2 ) terms cancel each other because of the 
specific coefficients in (|22p and the assumption (|81|) . Now, let us verify our proposal by explicitly 
calculating O (r) expansions of / (r, 7) and g (r, 7) . By using (f85|) in / (r, 7) , one obtains 



\ 



det 



(1 +7a ) 



\ 



det 



1 



+ 



7T 



ai 6 



1 + O (r 2 



, (87) 



1 J (1 +700) y ci di 
and it is possible to make the O (r) expansion by using (f22j) ; 

' ^ (ai + di) + O (r 2 ) =(l+7a ) + i7r(a 1 + d 1 )+o(r 2 ). 

(88) 



/(r, 7 ) = (I + 700) 



1 + 



2(l+7«o) 



Therefore, the O (7 2+n ) , n > 1 terms in 7 expansion of / (r, 7) do not contribute to the O (r) 
terms, only if (jSTjl holds. On the other hand, the O (r) expansion of o (r, 7) can be simply found 

as 



1 1 
9 ( T , 7) = 1 + ifl [2ao + r (ai + di)] = (1 + 7«o) + ipT («! + di) • 



(89) 



As a result, if one wants to consider O (r) behavior of / (r, 7), then one can equally work with just 
g (r, 7) which is simply equal to the O (7) expansion of / (r, 7) . 
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Appendix B: Analyzing Einstein-Hilbert Action and Quadratic Curvature Gravity with 

Second Order Perturbations 

In this Appendix, second order expansions of the curvature tensors are used in the well-known 
cases of the Einstein-Hilbert theory, and the quadratic actions including the Einstein-Gauss-Bonnet 
theory. This will help us construct the following O (/i 2 ) actions that frequently appear in the 
computations 

J d 4 x V^R(2), J d 4 x V=§gVR$, J d 4 x ^R^R^ U , 



J d*x ^g (Rl paX ) (2) , / d*x y/^Tt 1 { R W) . 



in terms of the building blocks appearing in Eq. (25) of 25]. 



(90) 



Analysis of the Einstein-Hilbert action 

First, let us find the second order in metric perturbation for Einstein-Hilbert action: 

I = - j d 4 x yf^j (R - 2A ) , 
and expanding up to third order in hnv yields 



(91) 



I=- <Tx W-g 

K 



1 + \h + V (h 2 - 2hl u ) + O (r 3 )| [(r - 2A ) + tR l + t 2 R {2) + O (r 3 ) 



- / d i xyj-g{ (R- 2A ) +r 



-h(R -2A Q ) +R L 



1 1 

R - 2A ) (h 2 - 2h%A + -hR L + R {2) 



+ T 



(92) 



One can find the nonlinear equation of motion for constant curvature background by investi- 
gating the first order term in r of the above action as 



~h(R-2\ {1 ) + R L 



after putting the explicit form of Rl and dropping out a boundary term one can get 



1 



0(h) 



drx y^gh(A - A ) , 



from which it follows that (A — Ao) g^v = upon taking variation with respect to huv- 
One can read the second order action as 



where R{2) can be read from (|13j) as 



o ( A ~ ) (Si^h - 2h pv ) + \g^R L 



+ R 



(2) 



R m = R px h ap ht - h^R%, - g ua K (R 



'.2. = n ni ,n x -ir n,„, - // w p yx' vtla ) L - g v °g m gfr 



(93) 



(94) 



(95) 



P ) 



(96) 



21 



Let us concentrate on f d 4 x y/—gR(2) part of the action and work out the integration by parts; 
J d 4 x V^R(2) = J d A x v^|a^ - h^R^ - t°h* (^V) L 



(n a ) L (Ti) L - {na)L (r 



(97) 



One can find g va h^ (R 13 ^) as 



4A 



(98) 



By using the definition of the linearized Christoffel connection in ([S]) and doing integration by 
parts, the last two terms in f d 4 x y/—g~R(2) can be found as 



d x sj — g 



(99) 



<i 4 x a/— 5 



(3DV 



Finally, J d 4 x \f—gR(2) becomes 

I d 4 x yf-lR(2) = I d 4 x yf = gh ilv (-^-RjL ~ ^9^Rl + A/t Mi/ - -g^^ > 



and putting this result in (|93|) yields 



and since A = Aq is found from equations of motion for constant curvature background; 



z O(/j 2 ) 



7L 



(100) 
(101) 

(102) 

(103) 



Analysis of the quadratic action 

Now, let us consider the quadratic actions in the form 

1 



d x y/—g 



(R - 2A ) + aR 2 + PR 



(104) 
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and calculate the second order action in metric perturbations. Then, up to third order, the expan- 
sion of the action is 



I 



j d A x V^I | - (R - 2A ) + aR 2 + $R\ 



1, fl 



R-2A ) +aR 2 + /3R 



+ [-Rl + 2aRR L + (3R fn/ R^ + (3 



(105) 



h 2 ~ 2h%, 



R-2A ) + aR 2 + (3R 



1, (\ 



( -R L + 2aRR L + m^Rjlu + P (R 



+ I -R {2) + 2aRR {2) + aR 2 L + PBTR® + (3 (R^)(i) R»v + P (^)( 2 ) R. 



'(1) R^v 



Here, note that (R^ u )^ and (R^)^) are the first and the second order terms in the metric 
perturbation expansion of R^ u . First of all, in order to find the nonlinear equation of motion for 
constant curvature background, one needs to study the O (r) term in the above action. After using 
the definitions of R^ v , Rl and dropping out the boundary terms one can get 



O(h) 



(fx y^gh(A - A ) , 



(106) 



which yields the equation of motion A = Ao- Then, let us move to the second order term in metric 
perturbation. After using the result given in f 1 1 1 [) . one can obtain 



z O(/i 2 ) 



- 7, { [ - + 8aA + 4/3A j W'Q 



L 



-¥ 



i 



(A - A ) + 2/3A 2 



1 



(A - A ) + 6/3A s 



(107) 



+2ah^ (g^D - + Ag^) R L 

-2(3 (Ag^R$ + R%R%, + R^Ag^) } . 

Here, let us first handle f d 4 x yj —gR^ R^ v . Using the definition of R^ and using the linearized 
Bianchi identity (and also its covariant derivative) which is 



i 



7 iiv 0, 

one can find the following result after use of integration by parts 



(108) 



d x 



gR^R^ u 



d A x V^gh^ 



[g^U - V fl V u + Ag^j R L + UQ t 



2A 



9^Rl 



14A 



A 



8A 2 



-Rjiv + -^9i-wRl + 



(109) 



( 2 ) 



Second, (R^)^ is related to R^J in the following way: 



(Rn { 2) = fa (g m 9 vP R a p) m = gTR$ - ^hTR^ + 3Ah 



2 



(110) 
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- (2) 

The g^Rfni term can be given in terms of R^) with 



R { 2) = {g^R^) m = sTR$ - ^R% + AC- 



(111) 



Then, f d A x yf—gg^Rpu becomes 

J d* x yf=tgTR$ = hr {^Riu - \g^RL - jg^h 

with the help of fjlOlj) . By use of these results and the equation of motion for constant curvature 
background which is A = Aq in I ^ h 2^, one can get 



(112) 



l O(h 2 ) 



2 1 ** 



- + 8aA + ^/3A ) Q 



7 /IV 



+ (2a + 0) (g^U - V M V, + Kg^) R L + (3 ( ug^ u - ^g^ v R L 



(113) 



which is same as Eq. (25) of 25]. 

Finally, let us analyze the Einstein-Gauss-Bonnet theory, 



I 



d x \J—g 



-(R-2A ) +1 (Rl pcTX -4Rl„ + R : 



(114) 



just to check the consistency of our construction. Here, the only remaining part that we have not 

? 2 
ppoX 

fipaX ' 



analyzed is the R p \ pa x term. First, let us use the previous result in order to obtain the second order 
action in metric perturbations for the terms other than R 2 



- + ^ 7 A^ ^ - 2 7 (^D - V^V, + A^) i? L 



4 7 



2A 



-gp,vRh 



<"" 3 

Then, up to third order, expansion of the last term becomes 



(115) 



/ 



7 = 7 d^xV-9 {R'x + r 



E>2 

^ppaX 



(1) + 1 / ) R 2 



(!) . 1 



(116) 



First of all, it should be shown that first order part is a boundary term such that it should not 
give a contribution to equation of motion for constant curvature background: 



0(h) 



where 



Then, 



R 



d x y/—g 



8A 2 



(Rl P «x) (1) + \hR 



fipaX 



ppaX 



3 ' 



>« . 1, „2 
{R^paxf ) = ItRl- 



l O(h) 



d x \J — g 



;ii7) 



(118) 



;ii9) 
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and since the remaining part is a boundary term, no contribution comes to the constant curvature 
background equation of motion from the square of the Riemann tensor. Then, moving to the part 
that is second order in metric perturbation 



Io(h*) = 7 / d x y/-g 



( 9 \(2) 2A A 2 / , 

(Rl paX y+—hR L + — (h 2 -2h 



(120) 



where 



2 

ppaX 



(2) 



E>M DP crp \-y 
tt pa\ a p-yu9 9 



(2) 



--2R p „ Xr7 ( R p 



pcrX 



(2) 



p<j\ 



(i) 



BP 



(1) 



+ 2 



R\\ (*Va) (1) + R%<? w sfi 1 ) + R^xR'^gft, (121) 



Using R^pa = f {g^gua - g^g vp ) and i?( 2) = g pa Rp 2 } + g p ^R P a + i^ff 



(2)' 



(i) 4A 



2A 2 



and using (198 



2 

pi/ / 

122) 



(i) 



.a) 



4A T 



14A 2 



2A 2 



pv 



Now, let us consider J d 4 a; y 7 — g \ R 2 ^ pa \ 
be found as 



(2) 



. The jd^x^gg^g^ (R^x)^ (R"^' term can 



(123) 



(i? A ffp J (1) (iT w )« =^ 



2A 



2 ( - y^ R l ) + " V M V^ + A^J i? L 

-hP* (30J2^ - 9^i? x - 32A V + 2A^/t) , (124) 



after a somewhat lengthy calculation where the definition of the linearized Riemann tensor is used 
and the terms are rearranged by using integration by parts. Using this result with (|ll)ip . one get 



<Fx y/-g[R 



(2) 



<fx V^h^ 



2 { u Gpu ~ y^Rl ) "I - (.'/,„.□ - V /( V, + \ fipi, ) R L 



A 



80£, - ±R%, + 6A V - A^/j) }> , 



(125) 



and plugging it in the action: 



7 / (fx V^h^ 



8A 
~3~ 



Qpu + 2 (gpuO - V M V„ + A^) R L + 4 ( - 



(126) 

and considering this result with the part of the action coming from ^R 2 and —A^R 2 ^ terms given 
in (|115p . one finds that all the 7 terms vanish, and the Gauss-Bonnet term does not contribute to 
the equation of motion. 
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Appendix C: Linearization of the O (R 3 ) Action 

The following formulae are needed for the linearization of the O (ai?) 3 equations. The 
quadratic parts below already appeared in [25J, we reproduce them here for the sake of com- 
pleteness, the cubic parts are new. 

2A„ r A_ „ A 2 



5(aR^) AD V, 5(V fl V u R) = V fl V I/ R L , 5(aR) = aR L , 

S (rpr^R*) = 3A 2 Rf lu - 2A 3 V. 6 (R^R 2 a p) = 4A 2 i^ + 2A% V R L , (127) 
5 (R^R 2 ) = 16A a fl£, + 8A%„R L , 5 {RR P U R^ P ) = 8A 2 i?£„ + A%„R L - 4A 3 h^, 



1„ - 



4A„ r A 



4A 2 



5 (V^V u R a(3 ) = V^uRap - AV^V u h a(3 . 

Here, last two equations can related by using linearized Bianchi identity: 

1 



= 0, 



£V = R^u ~ ^uRl - Ah 



(128) 



Appendix D: Coefficients for the R— R Expansion 



Coefficients in the expansion of the square root of (J68J) are 



dR 



a 



(RM) 2(1 + aA)' 



(l + 4aA + 6a 2 A 2 + 4a 3 A : 



df_ 
dR 



0j 



{RM) 



2(1 + aA) z 



aA + 3a 2 A 2 + 3a 3 A 3 ) , 



d 2 f 



OR 2 



a 



2(1 + aA)' 



(l + 4aA + 6a 2 A' 



— rr (l + 4aA + 6a 2 A 2 + 4a 3 A 3 ) , 

4(l + aA) 6 V 1 

(129) 



d 2 f 



8R p a dR a R 



Q- 



2(1 + aA) 



a 2 W 



4(1 + aA) 



(l + aA) 2 <5^-a 2 A 2 <^ 
g (aA + 3a 2 A 2 + 3a 3 A 3 ) 2 , 
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df 
dRdR a p 



or 



(R,Rv ) 



2(1 + qA)' 



(ah + 3a 2 A 2 ) 



4(1 + aA) 

Coefficients in the expansion of the (1 73 [) are 



- a „ f 1 + 4aA + 6a 2 A 2 + 4a 3 A 3 ) f aA + 3a 2 A 2 + 3a 3 A 3 
1 + aA) 6 v M 



f(R,K) 



a 3 A 3 ( 1 



aA 



df 



dR 



{r,rz) 



2a 3 A 2 I 1 



3aA\ 



df 
dR% 



5a 3 A 2 p _ a 4 A 3 p 



S 2 / 

dR 2 



(j?,-R^) 



3a 3 . 9a 4 A 2 . . 

— A + ^, (130) 



df 



dRdR 



a 



3 A 



(rM) 



9a 4 A 2 
16 



9 2 / 



dR p a dR<$ 



a 3 A 
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